Decoherence induced by an ordered environment 
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To study the effect on the time evolution of a qubit, of an environment which undergoes a tran- 
sition to an ordered phase, we consider a qubit weakly coupled to a standard BCS superconductor. 
We find that ordering in the superconducting bath generically leads to an unfavorable faster-than- 
exponential decay of the qubit coherence. On the other hand, a coupling to the non-ordered sector of 
the bath can result in a drastic reduction of the qubit decoherence. Since these behaviors are endemic 
to the ordered phase, qubits can serve as useful probes of phase transitions in their environment. 



The past decade has seen tremendous activity devoted 
to developing experimentally viable qubits for quantum 
computing. These two-level systems are realized, for ex- 
ample, by directly using the charge or spin degrees of 
freedom of electrons and quantum dots [l[ or more com- 
plex entities like flux qubits [U, Q and Cooper boxes 0] ■ 
The utility of all these qubits for quantum computation 
is strongly limited by the influence of their environment 
which tends to destroy their quantum coherence. Conse- 
quently, a lot of recent theoretical studies have focused 
on ways and means of increasing the coherence time 
scales [5l4lO| . Another point of view on the interaction of 
small quantum systems with their surroundings consists 
in employing such 'qubits' as probes of the environment 
as typically done in Nuclear Magnetic Resonance spec- 
troscopy [Uj. This could be especially of use to investi- 
gate both equilibrium and non-equilibrium properties of 
cold atoms, for which standard thermodynamic measure- 
ments are difficult to carry out [l^] ■ 

The first theoretical studies of decoherence considered 
environments consisting of independent harmonic oscil- 
lators [l3| or spins [lj, 15 1. More recently, complex 
baths have been studied. In particular,intrabath inter- 
actions have been taken into account 0, @, EH E3 > raising 
the question of the possible influence of thermodynamic 
phases and transitions on the decoherence of the qubit. 
Numerous results have been obtained but the complete 
picture is far from clear. It has been shown that the 
mitigatin g im pact of intrabath interactions seen in many 
cases 3, l3| breaks down when the bath is in the vicin- 
ity of a continuous phase transition where the decoher- 



ence rate diverges on the disordered side [6|. Refsll6l. Il8l 
argued that symmetry brea king in the reservoir helps re- 
duce decoherence, while Rcf j 1 71 found a strong Gaussian 
decay of quantum coherence. However, all these works 
suffer from different drawbacks ranging from a complete 
neglect of low energy modes in Reflla to obtaining an 
order parameter independent behavior for the time evo- 
lution of the qubit in Refll7 . 

In this Letter, we revisit the problem of ordered baths 
to have a clearer understanding of their effect on qubits. 



We consider a superconducting bath described by the 
Bardeen Cooper Schrieffer (BCS) theory. Though the 
BCS Hamiltonian docs not capture the fluctuations in the 
disordered phase, it provides a very good description of 
the ordered phase. In line with Ref|fj|, where it was shown 
that the impact of the ordering on the qubit depends cru- 
cially on the relation between the qubit-bath interaction 
and the order parameter, we study different kinds of in- 
teractions between the qubit and the bath. The ordered 
superconducting phase is found to be characterized by a 
rich variety of behaviors not seen in the disordered phase, 
including faster-than-exponential decay of the coherence. 
The latter makes it unfavorable from the point of view 
of quantum computing. But, there are interesting excep- 
tional qubit states which decohere slower when the bath 
orders. The sensitivity of the qubit to the order in the 
bath, makes it a good probe of the transition in the bath. 

The combined system of the qubit and the supercon- 
ducting bath is described by the Hamiltonian 
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(r q -V + H B , 



(1) 



where cr q is the vector Pauli operator for the qubit, 
whose components are the usual 2x2 Pauli matrices, 
V is some bath vector operator that will be specified 
later, and Hb is the conventional BCS Hamiltonian 
Hb = J2ke Ek&k t &ke (Hi- The Bogoliubov operators 
a.k e are related to the electron annihilation and creation 
operators by 
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where c\. e creates an electron with momentum k and 
spin e =t, I- The BCS dispersion relation is E k = 
sgn{ek)\/ e\ + A 2 , where e k is the underlying electronic 
dispersion and A is the superconducting gap. The co- 
efficients in © obey {u,v)\ = (1 ± e k /E k )/2. We set 
H = Icb = 1 in the rest of the paper. The superconducting 
order parameter A at temperature T is self-consistcntly 
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determined by 

y io Ve 2 + A 2 V ; 

where g is the strength of the phonon-mediated electron- 
electron interaction, ujjj is the Dcbyc frequency and N is 
the electronic density of states at the Fermi surface. As 
is well known, this equation determines a critical tem- 
perature T c which separates a high-temperature metallic 
phase where A = and a low-temperature phase where 
A increases monotonically as T decreases. We assume 
that, at time t = 0, the qubit and the bath are uncorre- 
cted and that the bath is in thermal equilibrium. The 
initial state of the combined system is thus ft = p(0)®pB 
where ps oc exp(— Hb/T) and p(0) is any qubit density 
matrix. The time evolution of the reduced density ma- 
trix of the qubit is given by p(t) = Trsf2(f) where Tib 
denotes the partial trace over the bath degrees of free- 
dom and f2(t) = exp(—itH)rtexp(itH). In the limit of 
weak coupling between the qubit and the bath that we 
are interested in, one can use two main methods to cal- 
culate p(t): i) the time-convolutionlcss (TCL) projection 
operator technique and ii) the Nakajima-Zwanzig (NZ) 
approximation [iol |2lj |. TCL gives local- in-time equa- 
tions of motion for pit) whereas the NZ approximation 
gives an integro-diffcrcntial dynamical equation for pit). 
The accuracy of these methods depends on the problem 
studied and it is difficult to assert a priori which one is 
more appropriate [22j |. In this Letter, we focus on the 
asymptotic evolution predicted by the second-order TCL 
approximation, and briefly discuss the results obtained 
using the NZ technique at the end. To write the the mas- 
ter equation given by the TCL approximation to second 
order, it is convenient to first rewrite the Hamiltonian ([T]) 
as H = H B + (T q ■ (V) + Hi where (...)= Tr(p B . . .) and 
Hi = er q ■ (V — (V)). With these notations, we obtain 

d t p = -i[a q -(V),p(t)}- f d-rTr B [Hi, [H t (-t), p(t)®p B ]], 
Jo 

(4) 

where Hj(t) = cxp(itHs)Hj cxp(— iti^g) = tr q ■ (V(i) — 
(V)). This equation is our starting point for the follow- 
ing calculations. We remark that, replacing the integral 
boundary t in ((H) by +oo leads to the well-known Marko- 
vian master equation. However, as we will show below, 
the time evolution of pit) can be non-Markovian in the 
ordered phase precluding the use of such Markovian mas- 
ter equations. 

Kondo coupling - We first consider a Kondo like cou- 
pling where the qubit couples to the electronic spin 
density at the origin S(0), i.e., V a = A5 a (0) = 
^J2k,k',e,e' c ke a ?e' c k'e' where a g {x,y, z}, A is the cou- 
pling strength, and o~f e , are the matrix elements of the 
Pauli matrix a a . The Hamiltonian (flj is effectively the 
same as that of a magnetic impurity embedded in a su- 
perconductor. Contrary to the situation of a metallic 



bath where one really does not have a true weak coupling 
regime because of the dynamical Kondo effect, here we 
have a weak coupling regime Q ■ The isotropy of the to- 
tal Hamiltonian (fTJ) and the absence of any net moment 
in the bath lead to the following simplifications: since 
(V) = 0, the first term in (jlj) vanishes and all components 
of the effective spin-1/2 corresponding to the qubit, i.e., 
s a (t) = Tr(/9(t)i7q ), satisfy the same equation of motion. 
Consequently, the qubit evolution is characterized by a 
unique time function M defined by s a (t) = M(i)s a (0) 
0]- We see that with a Kondo like coupling, both deco- 
herence and relaxation exhibit the same time evolution. 
By writing the electron operators Cfe c in terms of the Bo- 
goliubov quasi-particle operators we find 

lnM(t) ~ -8A 2 J rfi / n( ^ 2 /2) V n, (5) 

with the dynamical spin structure factor S + given by 

S^iui) — J c?e/ ± (e, u>)p(e)p(e — u})n(e)n(u} — e), (6) 

where ^(e.w) = I±A 2 /e(e-w), p(e) = |e|(e 2 - 
A 2 ) -1 / 2 is the superconducting density of states, and 
n(e) is the Fermi function. 

For T > T c , the bath is a simple metal (/ ± (e) = f) 
and we obtain the usual asymptotic Markovian decay 
lnM(i) = — jt with a decay rate 7 = 47rA 2 S' + (0) which 
increases and saturates to a density of states dependent 
value at very high temperatures @, H| . As one approaches 
the transition temperature T — > T c + , we expect the grow- 
ing fluctuations to result in a divergent rate 7 at T c Q 
though this is not captured by the mean field BCS the- 
ory used here. We now analyze the asymptotic qubit 
evolution in the ordered phase 0<T<X' c . At T = 0, 
we find that S + ius) = for all u> > —2 A, where 2A is 
the gap to two particle excitations. This gap leads to an 
incomplete (power-law) decoherence of the qubit just as 
seen in the case of an insulating bath 0. At tempera- 
ture < T < T c , due to the divergence present in the 
superconducting density of states, S + is infra-red diver- 
gent: S + (uj) ~ -r(T) In \u/T\ as uj -> 0. The pre-factor 
r(T) = A/2 cosh 2 (A/2T), shown in the inset of FigUJ is 
a non-monotonic function of temperature which vanishes 
at T = and T = T c . This infrared divergence results in 

In M(t) ~ -27rA 2 r (T)t In t (7) 

for times t > t a = Max(f/T, f/A). We see that, con- 
trary to naive expectations, the ordered bath leads to 
a novel faster-than-exponential loss of coherence of the 
qubit. This can be attributed to the fact that the qubit 
couples to the spin fluctuations and hence order param- 
eter fluctuations via the singlet Cooper pairs. Moreover, 
we see an interesting recntrance in the asymptotic regime 
because the coefficient r(T) which dictates the asymp- 
totic decoherence is the same for two different tempera- 
tures (cf Fig. [T]) . 
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FIG. 1: In M(t)/t as a function of lnt for T = 0.02, 0.03, 0.04, 
0.045 and 0.051 in units of uj). The inset shows the coefficient r as 
a function of T. Here gN = 0.33, T c ~ 0.056cj o and A(0) ~ O.lwrj. 



Order coupling - We now consider the coupling dc- 



0. 



v x -- 



V+V^ and V v = iV - iV^ 



scribed by V z = 

where V = AJ^ fe k , CjL.cL, is the order parameter at the 
origin. Such couplings are relevant for qubits made from 
Cooper pair boxes pi In this case, we obtain, from 
dU) , the following set of coupled differential equations for 
the components s a : 
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where E±(t) = —8\ 2 fdu>sm.(u)t)S ± (u))/u) and h = 
—2 (V x ) = 2AA/ gN. Unlike the Kondo case previously 
studied, here the equations for the components y and z 
are coupled in the ordered phase where A^0. A closer 
look shows that XL and hence s x are related to the dy- 
namic charge correlation function which is non-singular 
and £+ to the spin correlation function which shows sin- 
gular behavior. As will be discussed below, this leads to 
a complex time evolution of the qubit. 

For T > T c , all the components are uncoupled and 
decay asymptotically as lns Q = — ^ a t with the rates 
j x = -fy = 8it\ 2 S~(0) and ^ z = 27^. The existence 
of two different rates is a direct consequence of the spin 
anisotropy of the order coupling whereas in the spin 
isotropic Kondo case, all rates are the same. In the or- 
dered phase, at T = 0, solving Eq0 we find that the 
presence of a gap, both in E_ and S + , leads to an in- 
complete decay of the central spin. For T ^ 0, S~ is 
regular at low frequencies which results in a Markovian 
decay lns x ~ —~f x t for times t 3> t a . The full tempera- 
ture dependence of the rate ^ x is shown in the inset of 
Fig. H Its behavior for T -> is j x cx T e - A / T . This 
is an important result of our work because it shows that 
in the ordered phase, the Markovian rate for the com- 
ponent s x is strongly suppressed compared to a simple 
metallic bath. This is very similar to the relaxation in- 
duced by a coupling to the charge fluctuations, studied 
in the context of NMR by Fuldc and Black [25|, which 



is described by cr q ■ V cx a z h where n is the number of 
electrons at the origin. For T >T C , the rate j x coincides 
with the rate of the equivalent metallic bath with A = 
and saturates to a finite DOS-dependent value propor- 
tional to J dEp(E) 2 as expected. If the qubit's entire 
evolution was determined by this component, then the 
bath is effectively a semiconductor with a temperature 
dependent gap [8|. On the other hand, we find that the 
components s y and s z exhibit non-Markovian behavior. 
In the limit of weak coupling between th e q ubit and the 



bath, Eq.© can be solved for s y and s z 



By writing 



(s y ,s z ) = (s y (0), s z (0)) exp(A) and keeping only order 
terms of the 2x2 matrix A, we find 

s y ~e ai s y (0) (ansincO + cos 6) + s z (0)hisinc6 (9) 
- s v (0)ht sine 6 + s z (0) (cos 6 — aosincO) 



where a^t) = f* dt'\pT, + (t') + (ji - l/2)E_(t')], ©(t) 



211/2 



[{htf - a (t) 
that both s v and s 



ht, and sincO = _1 sinO. Note 
show oscillatory behavior in the or- 
dered phase with a frequency proportional to the order 
parameter. At T = 0, s y and s z oscillate at a maxi- 
mal frequency and at finite temperatures < T < T c , 
these oscillations are damped with the envelope a^it) cx 
—r(T)th\t for t ^> t a and we recover the faster-than- 
cxponential decay (|7|) with reentrant temperatures en- 
countered for the Kondo coupling (cf Fig. [J). The ap- 
pearance of oscillations of the qubit can be used as a tool 
to demarcate the phase diagram, measure the gap and 
the critical temperature of the superconductor. 

Because of the very different asymptotic behaviors of 
the components s a found above, the decay timescale of 
the qubit depends crucially on the initial conditions. As- 
sume the qubit is initially prepared in an cigenstate of 
a x . Then the components s y (t) = s z (t) = are con- 
stant, and the asymptotic time evolution of the qubit, 
determined by s x (t), is non-oscillatory and Markovian 
with a highly reduced rate in the ordered phase. These 
pure states are thus relatively stable in the environment 
considered here. Moreover, a combination of good initial 
preparation and pulse sequences which repeatedly orient 
the qubit in the x direction could efficiently reduce deco- 



herence times |27|,[28||. For an initial qubit state such that 
s x (0) = 0, this component remains zero and the asymp- 
totic behavior of the qubit consists of an oscillatory and 
faster-than-exponential decay. In this case, due to the 
coupling between y and z components, we have situations 
where even if one of the components s y or s z is initially 
zero, this component can grow with time and show oscil- 
latory behavior. For generic initial states, since s y and s z 
decay much faster than s x , the asymptotic time evolution 
of the qubit is Markovian. The qubit state first decoheres 
into a statistical mixture of the eigenstates of a x and then 
relaxes exponentially into the maximally mixed state. 
In both the Kondo and order coupling cases, we find 
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FIG. 2: s a (t) as a function of t for the order coupling to a super- 
conducting bath at T = 0.045w£> for the initial conditions indicated 
in the legend. For the components s y and s z , only the envelope 
is plotted. The inset shows 7^ as a function of T for metallic and 
superconducting baths. Here gN = 0.33 and A = 0.05. 



interesting intermediate time behaviors but these will be 
discussed elsewhere {29j. We now briefly discuss the re- 
sults obtained using the NZ approximation 0] . The con- 
clusions for the Markovian behavior seen for the com- 
ponent s x in the order coupling case remain unchanged. 
However, the asymptotic faster-than-exponential behav- 
ior seen in both the Kondo and order coupling cases 
is replaced by a much slower non-Markovian behav- 
ior Mjvz(t) ~ / dwcos(wt)/ln|a;| ~ 1/tlnt for times 
t 3> t nz . We find t nz 3> t a and an intermediate regime 
t a <C t <C t n z characterized by a quantitatively faster de- 
cay than that predicted by the TCL approach, where the 
qubit becomes practically incoherent. Anomalously fast 
decoherence seems to be a feature of both NZ and TCL 
methods in the ordered phase. 

To summarize, we have studied the influence of a true 
long-range ordered bath on the state of a qubit. We 
found, for two different qubit-bath couplings, a faster- 
than-exponential decoherence of the qubit in the ordered 
phase leading to the conclusion that ordered baths are 
disastrous for qubits. However, for one of the couplings 
considered, some particular pure states of the qubit are 
better preserved from decoherence when the bath orders, 
making such states potentially useful in experiments us- 
ing echo sequences or other related techniques. This re- 
sult constitutes another example of the important role 
played by the form of the coupling to a bath that can 
order Q. We expect, if fluctuations in the disordered 
phase are properly taken into account, a divergence of 
the Markovian rate at the transition from the disordered 
side, followed by a faster-than-exponential loss of coher- 
ence in the ordered phase for generic initial states of the 
qubit. Exceptions are possible, as our study shows. We 
believe that this picture should be valid for all ordered 
baths provided the qubit couples in some way to the or- 
der parameter. It would be interesting to generalize our 
model to a non-degenerate qubit and to study the effect 
of such anomalous dissipation on the tunneling of the 
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